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^~~^ • Abstract 

>; 
Cn I For different models of the electron-phonon interaction, the asymptotic 

^— V ' behaviour of the moments of the stationary homogeneous solution of the linear 

^SJ . Boltzmann equation is determined in the limit of a high external field. For 

^p I Hilbert-Schmidt kernels of a finite rank, a result recently proven for kernels 

Qv ' of rank one is found generally valid; as a consequence velocity saturation is 

"^ ■ excluded for these collision models. For a class of singular collision kernels in 

2 . contrast, velocity saturation is generally obtained. 

S; 

^b ■ 
d ; 1 Introduction 

o. 
o ■ 

IL" ■ It is a well known fact that the drift velocity of electrons in a semiconductor ap- 

• i— I . proaches a finite value, if the applied field becomes very large. Already Shockley, 

rS I whose experiments with RYDER [|I| first revealed saturation of the electron drift ve- 

cd I locity in germanium, argued in an intuitive manner that optical phonon scattering 

should be responsible for this effect [^ . Meanwhile, changeover from acoustic to op- 
tic mode phonon emission in high electric fields has been examined experimentally 
in some detail [@. 

For a mathematical explanation of velocity saturation, it is necessary to deter- 
mine the high field behaviour of the first moment with respect to the momentum p 
of the solution f{p) of the stationary homogeneous Boltzmann equation |^ 

(1) ^^^(P) = 2^(P) 

For several specific models of the electron-phonon collision operator Q, this be- 
haviour has been determined by different authors ||^, 0, 0. Particularly Reik and 
RiSKEN have shown that in the presence of both optic and acoustic mode deforma- 
tion potential scattering the electron drift velocity saturates H. Nevertheless the 



interesting question remains as to which general classes of collision models actually 
do or do not lead to velocity saturation. 

In this article we answer this question for two classes of collision models within 
the framework of a parabolic band approximation and the assumption of a nonde- 
generate electron gas, which means that the electron concentration be not too high 
IP). The main results are given in form of theorems at the end of section three and 
five. 

In section three and four, we examine regular kernels: if their rank is finite, 
the high field behaviour of the moments is given by a law which the author has 
proven recently for "relaxation time approximations" (kernels of rank one) . As 
a consequence, velocity saturation is excluded for these kernels. This result seems 
to remain valid for kernels of infinite rank as well, as we point out in section three. 

In section five we shall see that the situation changes considerably, if the quantum 
nature of the lattice vibrations is taken into account. Under the assumption of a 
single phonon frequency ("Einstein's model", which is a good approximation for 
the optical phonon mode [0) and some isotropy hypothesis upon the scattering 
probability, it will turn out that the electron drift velocity saturates, regardless of 
the detailed functional form of the scattering probability. 

2 The collision operator 

If the electron concentration is low enough, so that interactions among the electrons 
themselves (via collisions or Pauli's exclusion principle) may be neglected, the 
electron phonon collision operator Q in (1) is linear and has the general form |^ 
(throughout this article we set the effective mass of the electrons equal unity, so 
that momentum p and velocity are the same) 

(2) Qf{p) = M{p) J dp' K{p\p) f{p') - f{p) J dp' M{p') K{p,p') 

where M{p) := (27r0)^^/^e~^ /^^ is the Maxwellian at temperature 9 of the semicon- 
ductor lattice and the collision kernel K{p,p') is a symmetric and positive distri- 
bution. The second integral in (2), which gives the total scattering rate, is usually 
called collision frequency v 

(3) v{p) := / dp M{j)') K{p, p) 

The collision kernel K{p,p') strongly depends on the model for the electron-lattice 
interaction. We can roughly divide these models into two classes: 

The first class are classical models, which allow an arbitrary amount of energy 
exchange between electrons and lattice; these models lead to continuous collision 
kernels or (under the not very restrictive assumption that K be squareintegrable 



with weight M) Hilbert- Schmidt kernels. Nearly all existence and uniqueness results 
for the linear Boltzmann equation refer to collision kernels belonging to this class. 

The second class are semiclassical models, which take into account the quantum 
nature of the lattice vibrations. These models lead to singular kernels, where energy 
conservation is represented by Dirac (5-distributions in the collision kernel. 

3 Hilbert-Schmidt kernels of finite rank 

Let us begin our examinations with collision kernels that are measurable functions 
with a finite norm 

(4) \\K\\Ij := JdpJdp'M{p) M{p') K\p,p') < oo 



Since K is symmetric, the kernel K'{p,p') = J M{p)K{p,p')J M{p') is symmetric 
and squareintegrable. Hence it follows from standard Hilbert-Schmidt theory |1^ 



that K' has an eigenfunction representation which converges in L^(1R ). If we write 
this representation in terms of K, we find 

oo 

(5) K{p,p) = Y.^n ^n{p) ^n{p) 

n=0 

where the right hand side converges with respect to the norm (4). The n-th eigen- 
function (fn is a solution of the eigenvalue equation with eigenvalue A„ 

(6) J dp M{p')K{p,p)ipnip') = Xn^nip) 

The eigenfunctions {(pn} form an orthonormal system with respect to the L^-scalar- 
product with weight M: 

(7) dpifinip) Vmip) M{p) = 6nm 

Now let us assume that the eigenvalue equation (6) has a solution only for a finite 
number of eigenvalues {Xn}n=T y which means that the kernel K{p,p') is of finite 
rank. In this case the eigenfunction representation (5) takes the form 

(8) K{p,p)= J2 XnVn{,p)Vn{.p) 

n=l 

For the special case Umax = 1? (8) is identical with the relaxation time approximation 
discussed in |^. We shall see that the stationary homogeneous Boltzmann equation 
can be readily integrated for any finite rimax and the resulting expression for the 
moments allows an asymptotic evaluation for large fields E. 



3.1 Solution of the Boltzmann equation 

With the coUision kernel (8), the homogeneous stationary Boltzmann equation takes 
the form 

(9) [Edp+ J2 l^nVn{p))f\p) = M{p) J2 anVnip) 

n=l n=l 

where we have introduced the abbreviations 

oo oo 

(10) Hn ■= A„ / dp M{p)ifn{p) and a„ := A„ / dp f{p)ipn{p) 

— oo — oo 

From its definition (3) we see that the collision frequency in this model is given by 
^(p) = J2 fJ'n '^n{p) ■ The solution of (9) vanishing at p = ±oo is 



1 f^max p 

(11) fip) = -i e-^(^)/^ E a„ / dp' vn{p') M{p') e^(^')/^ 

"=1 -oo 

N{p) (read "capital z/") denotes the integral over the collision frequency 

(12) whereN{p) := / dqu{q) = dq Y] ^in'~Pn{q) 

Jo Jo „=i 

The coefficients a„ must be determined selfconsistently by insertion of the solution 
(11) into the definition of the coefficients (10). This leads to the system of linear 
equations 

ra=l 

where 6mn denotes Kronecker's symbol and the matrix A^n is given by 

oo p 

A™„ := - |rfpe-^(^)/^A„v^„(p) J dp' Mp') M{p') e^^^P'^^ 

— oo — oo 

Partial integration in the p-integral using N'{p) = J2 f^k'^kip) reveals that 

Yl 'V i-^kn - 4n ) = 

k=i Afc 

so that the determinant of the matrix (A — 1) vanishes. Hence the homogeneous 
system (13) indeed has a nontrivial solution and the stationary homogeneous Boltz- 
mann equation is solved by (11). 



3.2 Asymptotic evaluation of the moments 

After changing the order of integrations we obtain for the n-th moment of / given 
by (11) (remember that the solutions a^ of (13) depend on E) 



4"(E) := j dpe-^''^'^k{p)e''^^^''' jdqq'^e 



p 



If we make the substitutions t = N{q) and s = N{p), we can interpret the double 
integrals in (14) as Laplace transforms with the Laplace variable e := 1/E (we 
denote the inverse function of A^ by A^^^): 

(15) = Jdte- /* ^'='f;'-» '^-;'^+ -»'• e-(-'"»V- 

^ ^ J J z/(A^-i(s)) z/(A^-i(t + s)) 

According to the Abelian theorems on Laplace transforms |TI|, the asymptotic be- 
haviour of the Laplace transform (15) for small e is determined by the asymptotic be- 
haviour of the inner integral for large t. Under the assumption that N~^{t) behaves 
asymptotically as given in (16) (see theorem 1 below), the asymptotic behaviour 
of the s-integral for large t can be determined in complete analogy to section 4 of 
reference 0]. The result is 

s-mtegral m (15) ~ ^ttt-^ttt / ds — ,;, -,, ,, e ^ ^ '^ ' 
^ ' v{N-^{t))Loo z/(A^-i(s)) 

For small values of e = l/i?, we may replace the s-integral in (15) by this asymptotic 
expression for large t. If we resubstitute t = N{p) and insert the result into (14), 
we obtain 



oo , , , oo 






E^oo " 



Efc=r ak{E) f- jdpe-^(-py^ jdpe-^(py^ 

'' 



This is exactly the same expression that was obtained in for kernels of rank one in 
the zero temperature limit. Consequently we obtain the same high field behaviour 
of the moments (for details of the derivation see [0): 



Theorem 1: The collision kernel be of the form (8) and the inverse function N~'^ 
of the integral over the collision frequency N{p) = Jq v{q) dq have the asymptotic 
behaviour 

(16) N-\t)^ Af'Lit) for t ^ oc 

where a > —1 and L be a continuous function with the property limt^oo L{Tt)/L{t) = 
1 for every r > 0. 

Then the n-th moment of the solution of the stationary homogeneous Boltzmann 
equation behaves for E ^ oo like 

(17) {p-) ~ ^V^";; + ^ + ^^ i^-L"(E) ~ const. {N~\E)y 

n + 1 T{a + 1) ^ ^ 

In particular the first moment or the drift velocity behaves as N~^{E) for large 
E. Because of N{oo) = cxd, which is a necessary condition for the existence of a 
stationary homogeneous solution of the Boltzmann equation |jl2|, we conclude that 



velocity saturation is excluded for Hilbert Schmidt kernels of finite rank. 

4 A conjecture 

Since we have seen that theorem 1 on the high field behaviour of the moments is 
valid for Hilbert-Schmidt kernels of any finite rank, it is reasonable to assume that 
our theorem holds for arbitrary Hilbert-Schmidt kernels as well. This conjecture 
however is very difficult to prove; a major problem is that the eigenfunctions are 
not positive for a rank greater than one, so that their sum (5) is not necessarily 
positive if it is truncated at finite n. Maybe a proof is possible, if the kernel K{p,p') 
is a continuous function of p and p', for then the sum (5) converges even uniformly. 
Nevertheless the problem mentioned above still remains. 

As a reinforcement of our conjecture, let us consider the kernel of infinite rank 
K{p,p') = p\p — p'\, for which Gervois and PlASECKl determined the large field 
behaviour of the first moment in B . They obtained 



(18) {p){E)^ ^2E/p7i ioTE^oo 

which is independent of the temperature 6. In this model, the collision frequency 
reads 



u{p) = 2pe M{p) + p\p\ il-2 r M{p') dp 



'\p\ 
Hence, the integral over the collision frequency behaves as N{p) ~ pp'^/2 for large p 



and its inverse function as A^ ^(^) ~ v 2t/p for large t. Naive application of theorem 
1 leads to 

{p)iE) -X- v^ 2r(0.5^'l) ^ = ^/^^^ 

6 



which is identical with Gervois' and PlASECKl's result (18)! 

5 Optical phonon scattering 

Now let us consider collision models that take into account the quantum nature 
of the lattice vibrations. Under the assumption of a single phonon energy uj (this 
assumption holds for the optical mode), the general form of the collision kernel reads 



(19) M{p)K{p\p) = (3{p\p){{No + l)6{p'/2-p'^/2 + uo) 

+ No5{p'/2-p'^/2-uj)] 

where No = {e^^^ — 1)^^ is the occupation number of the optical phonons with energy 
u and the scattering cross section f3{p',p) is a symmetric and positive function. 
The first term in the curly braces describes phonon emission, the second phonon 
absorption during the scattering process; the Dirac (5-distributions assure energy 
conservation. 

Under the isotropy assumption /? = /?(|p'|, \p\), we shall see that an asymptotic 
solution of the Boltzmann equation is possible for high energy electrons with p'^^2uj, 
provided the electric field is sufficiently high E^uj. The asymptotic knowledge of 
the electron distribution will be sufficient for the determination of the high field 
behaviour of the moments. 

Our replacement of difference operators by differential operators, which is the 
essential step in our procedure, is not new. The first who did so were Yamashita 
and Watanabe H^. However, they obtained erroneous results for the electron 



distribution and the drift velocity, which were corrected by Reik and RiSKEN P]. 
Hansch and SCHWERIN have applied this method to inhomogeneous situations as 



well [15|. In all three references, a combination of acoustic and optical phonon 



scattering is considered and the scattering cross section j3 is taken to be constant. 

5.1 The stationary homogeneous Boltzmann equation 

Under the isotropy assumption [3{p\p) = P{\p'\, \p\), the one dimensional stationary 
homogeneous Boltzmann equation with the collision kernel (19) reads 

{Edp + {No + l)e(p2 _ 2u;)/3-{p) + Nof3+{p)}f{p) 
(20) = {No + 1) (3+{p) ^ |/ (-v/p2 + 2a;) + / ^+^p^ + 2a;) | 

+ Noe{p' - 2uj) f3- {p) \ {/ {-^|7^^) + / (+v/p'-2c.)} 



where G denotes Heaviside's step function, and the abbreviation f3^ means 



±2uj 



Now let us spht f{p) = fo{p) + feip) iiito its odd and even part foip) '■= Wip) ~ 
f{—p))/2 and fe{p) '■= {f{p) + /(^p))/2. Then we may write (20) as two coupled 
equations for fo and /e for p > 0. With the substitution q := p'^/2 and the notations 



(21) 



(22) 



Fo{q) := fo ( V2g ) = fo{p) = -[f{p) T f{-p) 



and 



B{q) := 13 J2{q - cu), J2q J q{q - u 



we obtain the coupled system of equations for Fo (remember that /9^ are even 
functions of p and that f3{p',p) is symmetric with respect to exchange of p and p') 

(23a) EdgFoiq) = {No + l)B{q + uj)F,{q + uo) + No^iq- uo)B{q)F,{q - uj) 

- \No + l)e(g - uo)B{q) + N,B{q + uo)\F,{q) 



(23b) Ed,F,{q) 



{No + l)0(g - uj)B{q) + NoB{q + u) Fo{q) 



Using (21), we may express the moments of the electron distribution in terms of 
the functions Fo 



(24) 



{p-){E) 



f-^dpp-f 
I - oo dp f 



for odd n 



rdq{2q)^F4q) 

°r^, n.^ N/ /o- ^OT even n 



5.2 Asymptotic solution for fast electrons 

For sufficiently high electric fields, an asymptotic solution of (23) in the region p^ ^ 
2uj is possible via replacement of the difference operators by differential operators. 
For this purpose we switch to the new varible x:=q/E and denote the resulting new 
functions with hats: 

Fo{x) := Fo{Ex) = Fo{q) and B{x) := B{Ex) = B{q) 
After a rearrangement of terms, the system (23) reads for x > e := uo/ E 
(25a) d^Fo{x) = {N^ + 1){b{x + e)Fe{x + e) - B{x)F,{x)} 

- No{b{x + e)Fe(x) - B{x)Fe{x - e)} 

(25b) dj,{x) = -{{No + l)B{x) + NoB{x + e)}Fo{x) 



For small e <^ 1 and large x ^ e (which means p^ ^ 2uj) , the difference operators 
on the right hand side can be replaced by differential operators. If we make a Taylor 
expansion of Fe{x ± e) and B{x + e) around x and keep only terms up to first order 
in e, we obtain 



(26a) d,Fo{x) ^ ed,[B{x)F,{x) 

(26b) dJeix) ^ -{{2No + l)B{x)+eNod,B{x)}Foix) 

Taking into account that Fo{oo) = 0, we may integrate (26a) from x to infinity and 
insert the result into (26b). Again we keep only terms up to first order in e. 

(27a) Fo{x) ^ eB{x)Fe{x) 

(27b) d^F,{x) ^ -e{2N, + l)B\x)F,{x) 

Equation (27b) easily can be integrated with the boundary condition -Fe(oo) = 0. 
Thus we obtain an asymptotic solution of (25) for x ^ e under the assumption that 
e = uj/ E be small. We may also take into account the region where x ^ t does 
not hold by introduction of two functions Qo{t) with the property limt_,oo Qo(t) = 1. 
Then the factors go{x/e) describe transition to the asymptotic domain x ^ e and 
the solution globally valid in x reads 

(28a) F,{x) -' Cgeix/e)exp(^e{2No + l)£dyB\y)^ 

(28b) Foix) -' Cgoix/e)eB{x)exp(^e{2No + l)l^^dyB\y)^ 

with some arbitrary constant C > 0. 

5.3 High field behaviour of the moments 

If we insert the high field solution (28) into the expression for the moments (24), we 
find after the substitutions q = Ex and q' = Ey (remember that e = u/E) 



(29) (p")(E)''i"/„(E)//o(i?) with 



oo 



uj{2No+l)A{q)\ 



f Jdq{2q) — go (^j B{q) exp (^_^£i™^L2^j for odd n 

(30) In{E) = { oo^ u^i2No + l)Aiq) \ r 

J dq (2g) 2 g^ ( ^ j exp ( ^^ — ^2 ) for even n 

where we have introduced the abbreviation 

(31) Aiq):= f'dq'B^q') 

Jo 

As E goes to infinity, the integrals in the numerator and denominator of (29) diverge, 
whilst the integrals remain finite if the upper limit is replaced by a finite value. 

9 



Hence their behaviour for large fields is determined by the asymptotic behaviour of 
the integrands for large q and consequently our knowledge of the asymptotic form 
of the electron distribution for p^ ^ 2uj will suffice for the determination of the high 
field moments. 

For convenience, let us assume that A{oo) = oo. Then we can procede as in sec- 
tion 3 by writing the integrals in (30) as Laplace transforms, which we can evaluate 
with the Abelian theorem again. To this purpose we make the substitution t = A{q) 
and obtain for large values of E (the inverse function of A is denoted by A~^) 



(32) 



In{E) 






~ < 



exp f- 



uj(2No+l)t 
£2 



B^A-\t) 



n-l 
2 



exp f- 



Lu{2No+l)t 



for odd n 



for even n 



Under the assumption that B{q) asymptotically grows like a power for large values 
of q, the asymptotic large field behaviour of the integrals in (32) is easily evaluated 
with the aid of the Abelian theorem on Laplace transforms ([^ part I, chapter 13 
§1 theorem 7). The elementary calculation leads to 



Theorem 2: The collision operator be of the form (19) and the scattering cross 
section (3 satisfy the isotropy property f3{p',p) = l3{\p'\,\p\) and be such that for 
q —^ CO 



(33) 



B{q) = /3 (v^2(g-cu), y^) /^q{q - oo) 



a ■ q 



with some constants a > and a > —1/2. 

Then the n-th moment of the solution of the stationary homogeneous Boltzmann 
equation behaves for E -^ oo like 



(34) 



(pliE) 



~ < 



22aa; 



r 



£( 



n+l+2a 
4a+2 



^ V4a- 



4a+2 



2a + 1 



a'io {2No + 1) 



n+2a 
4q + 2 



£'2q+i for odd n 



ra+l 
4a+2 



r 



1 



4a+2 



2a + 1 



a^LU (2No + 1) 



_E2a + l 



for even n 



Particularly, for the first moment we find 



(35) 



{P){E) 



r 



g+l 
2a+l 



'(4a + 2)cj 



r i-LJ] V 2No + l 



which is a finite value. Consequently the drift velocity saturates if the electrons 
interact with optical phonons! 



10 



5.4 An example 



For the particular choice l3{p,p') = C\p\ ■ \p'\ proposed by Mahan [J16l, the function 
B{q) defined by (22) is simply a constant, so that the stationary homogeneous 
Boltzmann equation can be solved exactly in the zero temperature limit No = 



T7|| . The resulting drift velocity is 
(36) {pHE) . ^'"^ 



2uj + e^'^^Ti/2z erfc {^/zuj) 

where z is the positive real root of the transcendental equation (zE)"^ = 1 — e~^'^. For 
large E, this solution is approximately z{E) ^ u/E"^ so that the current approaches 



the finite value jsat = J^uj/ir as E goes to infinity. This result is in agreement with 
(35) for A^o = and a = 0. 

In this model the asymptotic formula for the drift velocity is easily improved by 
keeping terms up to second order in (25). If we do so, we find for the first moment 



(37) {p){E)r. 2^ 



V 7r(l + u;2/2E2) 

A comparison of our asymptotic formula (37) with the exact value (36) is given in 
the figure. The agreement for large values of the electric field E is clearly visible. 
(The figure is omitted in the cond_mat-preprint). 

6 Summary and conclusions 

In this article we have seen that the experimenatlly observed saturation of the 
electron drift velocity in semiconductors cannot be obtained with regular collision 
kernels of finite rank. This result seems to remain valid even for infinite rank, 
though we could not present a rigorous proof of this conjecture. An extension of 
our argumentation to three dimensions is straightforward and can be performed in 
complete analogy to the case of rank one . 

In his investigation of the null space of the Boltzmann collision operator 



Majorana has already observed that singular collision kernels have mathematical 
properties rather different from regular collision kernels. In this article we have seen 
that also velocity saturation is a feature restricted to singular kernels. 

Under the assumption of isotropic optical phonon scattering we were able to 
derive velocity saturation, regardless of the detailed form of the scattering cross 
section. However two interesting questions remain, which are worth further investi- 
gation. First, is there a generalization of our proof to three dimensions beyond the 
framework of a two term Legendre approximation? The latter approximation, used 
in |]14| and f^, is an effective one dimensional approximation, which is equivalent 



to our one dimensional model. Second, does our result remain unchanged if the 
scattering is anisotropic? 

11 
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